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TITLE OF THE INVENTION 

ELECTROMAGNETIC SURVEYING FOR HYDROCARBON RESERVOIRS 
BACKGROUND OF THE INVENTION 

The invention relates to seafloor electromagnetic surveying for oil and other 
5 hydrocarbon reserves. 

Determining the response of the sub-surface strata within the earth's crust to 
electromagnetic fields is a valuable tool in the field of geophysical research. The 
geological processes occurring in thermally, hydrothermally or magmatically active 
regions can be studied. In addition, electromagnetic sounding techniques can provide 
10 valuable insights into the nature, and particularly the likely hydrocarbon content, of 
subterranean reservoirs in the context of subterranean oil exploration and surveying. 

Seismic techniques are often used during oil-exploration expeditions to 
identify the existence, location and extent of reservoirs in subterranean rock strata. 
Whilst seismic surveying is able to identify such structures, the technique is often 
1 5 unable to distinguish between the different possible compositions of pore fluids within 
them, especially for pore fluids which have similar mechanical properties. In the field 
of oil exploration, it is necessary to determine whether a previously identified 
reservoir contains oil or just aqueous pore fluids. To do this, an exploratory well is 
drilled to determine the contents of the reservoir. However, this is an expensive 
20 process, and one which provides no guarantee of reward. 

Whilst oil-filled and water-filled reservoirs are mechanically similar, they do 
possess significantly different electrical properties and these provide for the possibility 
of electromagnetic based discrimination testing. A known technique for 
electromagnetic probing of subterranean rock strata is the passive magneto-telluric 
25 (MT) method. The signal measured by a surface-based electromagnetic detector in 
response to electromagnetic (EM) fields generated naturally, such as within the earth's 
upper atmosphere, can provide details about the surrounding subterranean rock strata. 
However, for deep-sea surveys, all but those MT signals with periods corresponding 
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to several cycles per hour are screened from the seafloor by the highly conductive 
seawater. Whilst the long wavelength signals which do penetrate to the seafloor can 
be used for large scale undersea probing, they do not provide sufficient spatial 
resolution to examine the electrical properties of the typically relatively small scale 
5 subterranean reservoirs. Moreover, since MT surveying relies primarily on 
horizontally polarised EM fields, it is intrinsically insensitive to thin resistive layers. 

Nonetheless, measurements of electrical resistivity beneath the seafloor have 
traditionally played a crucial role in hydrocarbon exploration and reservoir assessment 
and development. In industry, sub-seafloor resistivity data have generally been 

10 obtained almost exclusively by wire-line logging of wells. There are, though, clear 
advantages to developing non-invasive geophysical methods capable of providing 
such information. Although inevitably such methods would be unable to provide 
comparable vertical resolution to wireline logging, the vast saving in terms of 
avoiding the costs of drilling test wells into structures that do not contain 

1 5 economically recoverable amounts of hydrocarbon would represent a major economic 
advantage. 

In research fields that are not of commercial interest, geophysical methods for 
mapping sub-seafloor resistivity variations by various forms of electromagnetic 
surveying have been under development for many years (e.g. Sinha et al., 1990; Evans 
20 et al., 1994). WO 00/13046 and WO 00/57555 make proposals for finding 
hydrocarbon reservoirs using such electromagnetic surveying. 
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SUMMARY OF THE INVENTION 

The invention discloses a new approach for electromagnetic surveying to 
locate hydrocarbon layers. New source-detector geometries are used based upon an 
5 electromagnetic source. 

According to a first aspect of the invention there is provided an 
electromagnetic survey method for surveying an area previously identified as 
potentially containing a subsea hydrocarbon reservoir, comprising: obtaining first and 
second survey data sets with an electromagnetic source aligned end-on relative to a 
10 first detector and aligned broadside relative to a second detector. 

The terms source and detector are used interchangeably with transmitter and 
receiver respectively throughout this document. 

The survey data from end-on and broadside alignments taken in combination 
allow the difference between galvanic and inductive contributions to the signals 
15 collected at the detector to be determined. Collecting survey data highlighting only the 
galvanic contribution in an end-on geometry is not reliable. As is demonstrated in the 
examples given below, it is generally impossible to differentiate between a rock 
formation containing a hydrocarbon reservoir and one which does not contain a 
hydrocarbon reservoir by studying the end-on survey data alone. The previously 
20 proposed electromagnetic survey methods for finding hydrocarbon reservoirs are thus 
believed to be at best highly unreliable. 

The end-on survey data are sensitive to the presence of resistive hydrocarbon 
layers (exploiting largely galvanic effects and vertical components of induced current 
flow). By contrast, the broadside survey data are sensitive to the larger scale structure, 
25 but relatively insensitive to resistive hydrocarbon layers (exploiting the dominantly 
inductive effects). The reason why collection of survey data from the broadside 
geometry is essential for reliable electromagnetic surveying is that many features other 
than hydrocarbon reservoirs can affect the resistivity beneath the seafloor and the 
results of a survey. For example, resistivity often increases steadily with depth in 
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submarine sedimentary basins, due to the progressive expulsion of pore fluids by 
rising overburden pressure. Such a resistivity profile will produce effectively the same 
response in the end-on survey data, as would presence of a hydrocarbon layer. 
Comparative examples given below demonstrate this effect. 
5 In the preferred implementation of the survey of the first aspect of the 

invention, the end-on and broadside alignments correspond to an arrangement of the 
electromagnetic source and the first and second detectors in which a right angle is 
formed between a first line leading from the first detector to the source and a second 
line leading from the second detector to the source, and wherein the source has its 

10 dipole axis aligned along the first line. However, in practice, an only approximate 
satisfaction of this condition will not greatly reduce the quality of the collected survey 
data. In any case, this ideal condition will not be satisfied exactly in practice, since the 
source is typically moved during surveying, being in the form of an antenna towed by 
an unmanned submarine craft. It would also be possible to obtain useable data if the 

15 above-mentioned right angle was changed to an angle away from 90 degrees, for 
example anywhere from 45-135 degrees may be satisfactory. How much the quality of 
the data deteriorates as the survey geometry moves away from right angles has not 
been studied, although this would be straightforward to do using the modelling 
techniques described herein. 

20 With the first aspect of the invention, the first and second survey data sets are 

preferably obtained concurrently. This can be achieved during a single tow of the 
electromagnetic source. The data collected by the detectors can then be time- 
synchronised to the same absolute clock. 

According to a second aspect of the invention there is provided an 

25 electromagnetic survey method for surveying an area previously identified as 
potentially containing a subsea hydrocarbon reservoir, comprising: obtaining first and 
second survey data sets with an electromagnetic source aligned end-on and broadside 
respectively relative to a first detector. 
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The method of the second aspect of the invention thus differs from that of the 
first aspect in that a single detector can be used to collect both the end-on and 
broadside survey data. This can be done by towing the source twice, once in a 
direction along the line connecting the source to the detector, and again in a direction 
5 transverse thereto. Namely, the first and second survey data sets can be obtained 
consecutively. The relative alignment between the source and detector when the end- 
on and broadside survey data are collected can be varied in the same way as discussed 
above in relation to the first aspect of the invention. In other words, it is best if the 
dipole is aligned along a line connecting the source and detector when the end-on 
10 survey data is collected, and aligned perpendicular to that line when the broadside 
survey data is collected. However, deviations from that condition will occur and can 
be tolerated. 

It will also be understood that the first survey data set can be obtained before 
or after second data survey set using the same electromagnetic source. 

15 In principle, the first and second survey data sets could be obtained with 

separate sources carried by different transmitters, in which case the data sets could be 
obtained concurrently. However, in practice, it is likely that only a single source will 
be deployed and the first and second data sets will thus be obtained one after the other. 
According to a third aspect of the invention there is provided a method of 

20 analysing results from an electromagnetic survey of an area potentially containing a 
subsea hydrocarbon reservoir, comprising: providing first and second survey data sets 
obtained from an electromagnetic source aligned respectively end-on and broadside 
relative to a detector; and combining the first and second survey data sets to obtain a 
results data set that represents a difference between the end-on and broadside 

25 alignments. 

The method can be greatly improved by normalising each of the first and 
second survey data sets relative to respective first and second normalisation data sets 
or functions specific to the end-on and broadside alignments respectively, prior to the 
combining. 
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The first and second normalisation data sets or functions can be calculated 
from a rock formation model, or from the first and second survey data sets. 

In a preferred embodiment, the first and second data sets each comprise radial 
and azimuthal components of electric field or magnetic field measured at the detector, 
5 and the method further comprises: transforming the radial and azimuthal components 
into at least one polarisation ellipse parameter, prior to the combining. The 
polarisation ellipse parameter(s) can be the amplitude and/or phase of the component 
of the electric field or magnetic field aligned along a major axis of the ellipse. 

The method may advantageously further comprise: visually representing the 
10 results data set in a plot of at least two dimensions corresponding to the survey area. 
The visual representation can be a two-dimensional (2D) plot in plan view, or a three- 
dimensional (3D) plot including depth, e.g. a perspective view of the survey area. 

The plot may include markings of areas of equal or similar electromagnetic 
field strength. These markings could be contour lines, or colour or grey scale 
15 gradations with one colour or grey tone being used for a range of electromagnetic field 
strength values. Stepwise colour gradation is used in the preferred embodiment, with 
each colour representing a defined range of data values. 

When the plot is of normalised survey data it is helpful if the plot also includes 
lines of equal absolute electromagnetic field strength, showing how the signal strength 
20 has decayed as one moves away from the source. The lines, which may appear as 
contours, can be labelled with either a relative or absolute decay value. 

According to a fourth aspect of the invention, there is provided a method of 
planning an electromagnetic survey of an area identified as potentially containing a 
subsea hydrocarbon reservoir, comprising: creating a model of the area to be surveyed, 
25 including a rock formation containing a hydrocarbon reservoir and a body of water 
above the rock formation; setting values for water depth, depth below the seafloor of 
the hydrocarbon reservoir, and resistivity structure of the rock formation; performing a 
simulation of an electromagnetic survey in the model of the survey area by calculating 
first and second survey data sets for an electromagnetic source aligned end-on and 
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broadside relative to a detector; and combining the first and second survey data sets to 
obtain a results data set that represents a difference between the end-on and broadside 
alignments. 

In typical use, the simulation will be repeated for a number of source 
5 frequencies and/or source-to-detector distances, and/or other parameters which can be 
varied during a survey. This iterative procedure can be used in order to select 
optimum surveying conditions in terms of source frequency and source-to-detector 
distance for probing the hydrocarbon reservoir. By optimum, it is not necessarily 
meant that the best survey conditions are found, but only that a set of survey 

10 conditions is found which will provide strong, unmistakable indications in the case 
that there is a hydrocarbon reservoir at the survey site. The iterative procedure may be 
purely under manual control. However, preferably, the simulator can allow the user 
the option of automatically optimising the survey conditions. The user can then switch 
between manual and automated iteration as desired. 

15 The model should preferably include a body of air above the body of water, so 

that the simulation can take account of signal propagation paths including the body of 
air when calculating the first and second survey data sets. The propagation path 
through the air (the 4 air wave') will in fact dominate for shallower water and longer 
distances between source and detector(s), as will be apparent from the examples 

20 described further below. For deep water and shorter detector-source distances the 
effect is less important and may be omitted from the model. 

The method preferably further comprises: normalising each of the first and 
second survey data sets relative to respective first and second normalisation data sets 
or functions specific to the end-on and broadside alignments respectively, prior to the 

25 combining. The first and second normalisation data sets or functions can be calculated 
from a rock formation model, for example. 

The method preferably further comprises visually representing the results data 
set in a plot of at least two dimensions corresponding to the survey area. The visual 
representation can be a 2D plot in plan view, or a 3D plot including depth, e.g. a 
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perspective view of the survey area. The other comments made above in relation to 
the plots of the third aspect of the invention also apply to the fourth aspect of the 
invention. 

Another aspect of the invention relates to a computer program product bearing 
5 machine readable instructions for implementing the method of analysing results from 
an electromagnetic survey according to the third aspect of the invention. 

Another aspect of the invention relates to a computer apparatus loaded with 
machine readable instructions for implementing the method of analysing results from 
an electromagnetic survey according to the third aspect of the invention. 
10 Another aspect of the invention relates to a computer program product bearing 

machine readable instructions for implementing the method of planning an 
electromagnetic survey according to the fourth aspect of the invention. 

Another aspect of the invention relates to a computer apparatus loaded with 
machine readable instructions for implementing the method of planning an 
1 5 electromagnetic survey according to the fourth aspect of the invention. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

For a better understanding of the invention and to show how the same may be 
carried into effect reference is now made by way of example to the accompanying 
5 drawings. 

Figure 1 shows a schematic representation of an exploratory EM sounding 

survey. 

Figure 2A is a graph showing the modelled electric field strength as a function 
of distance from, and along the axis of, the EM source. This is shown for several 
10 frequency components. The modelled sea depth is 750 m. 

Figure 2B is a graph showing the modelled electric field strength as a function 
of distance from, and along the axis of, the EM source. This is shown for several 
frequency components. The modelled sea depth is 900 m. 

Figure 2C is a graph showing the modelled electric field strength as a function 
15 of distance from, and along the axis of, the EM source. This is shown for several 
frequency components. The modelled sea depth is 1300 m. 

Figure 3 is a schematic plan view defining a survey geometry coordinate 

system. 

Figure 4A shows a schematic section view of a model subterranean strata 
20 configuration containing a hydrocarbon reservoir. 

Figure 4B shows a schematic section view of a model subterranean strata 
configuration containing layers of increasing resistivity. 

Figure 4C is a graph showing the electric field strength (E 0 , E90) as a function 
of distance from the source for two different detector positions for the subterranean 
25 strata shown in Figure 4A. Also shown are the normalising signals (N 0 , N90) expected 
from a uniform subterranean strata configuration. 

Figure 4D is a graph showing the electric field strength (Eo, E90) as a function 
of distance from the source for two different detector positions for the subterranean 
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strata shown in Figure 4B. Also shown are the normalising signals (No, N90) expected 
from a uniform subterranean strata configuration. 

Figure 4E is a graph showing the electric field strengths (Eo, E90) shown in 
Figure 4C normalised by the normalising field strengths (N 0 , N90) shown in the same. 
5 Figure 4F is a graph showing the electric field strengths (Eo, E90) shown in 

Figure 4D normalised by the normalising field strengths (N 0 , N90) shown in the same. 

Figure 5 is a schematic plan view showing an arrangement of sixteen detectors 
on a section of sea floor above a subterranean reservoir, a particular source transmitter 
tow path is also shown. 
10 Figure 6 A shows a schematic section view of a model subterranean strata 

configuration containing a hydrocarbon reservoir. 

Figure 6B is a plan view representing the normalised electric field strength as a 
function of location for the model subterranean strata configuration shown in Figure 
6A for a first transmitter position. 
15 Figure 6C is a plan view representing the normalised electric field strength as a 

function of location for the model subterranean strata configuration shown in Figure 
6A for a second transmitter position. 
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DETAILED DESCRIPTION 

Figure 1 of the accompanying drawing shows schematically a surface vessel 
14 undertaking EM sounding of subterranean rock strata 8 within which a 
hydrocarbon reservoir 12 is located. The surface vessel 14 floats on the surface 2 of 
the sea 4. A deep-towed vehicle 18 is attached to the surface vessel 14 by an umbilical 
cable 16 which provides an electrical, optical and mechanical connection between the 
deep-towed vehicle 18 and the surface vessel 14. The deep-towed vehicle 18 is towed 
by the surface vessel 14 such that it remains consistently close to the seafloor 6. This 
is facilitated by an echo-location package 20 which relays information about the 
height of the deep-towed vehicle 18 above the seafloor 6 to the surface vessel 14. The 
deep-towed vehicle 18 receives electrical power from the ship's on-board power 
supply via the umbilical cable 16. 

A cycloconverter unit 30 generates the chosen waveform to be supplied to an 
EM source in the form of a transmitting antenna 22 which is towed by the deep-towed 
vehicle 18. The transmitting antenna 22 broadcasts the EM signal into the sea 4, and 
this results in a component passing through the rock strata 8. One or more remote 
instrument packages 26 record the signal received by receiving antennae 24 in 
response to the transmitted EM signal. If the separation of the transmitting antenna 22 
and the receiving antenna 24 is greater than a few hundred meters, the highly 
conductive seawater strongly attenuates the direct signal between them. The 
components of the EM signal that have travelled through the rock strata 8 and the 
reservoir 12 dominate the received signal and provide information about the electrical 
properties of these regions. At the end of the sounding experiment, a remotely 
operable release system allows the instrument package 26 to be detached from a 
ballast weight (not shown) so that an in-built flotation device 28 can carry the 
instrument package 26 to the surface 2 for recovery and retrieval of data for inversion 
analysis. 
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The transmitting antenna 22 emits signals that propagate outwards both into 
the overlying water column 4 and downwards into the seafloor 6 and the underlying 
strata 8, 12. In both cases, at practical frequencies for this method and given the 
typical resistivity of the media 4, 8, 12, propagation occurs by diffusion of 
5 electromagnetic fields. The rate of decay in amplitude and the phase shift of the signal 
are controlled both by geometric spreading and by skin depth effects. Because in 
general the underlying strata 8, 12 are more resistive than seawater 4, skin depths in 
the underlying strata 8, 12 are longer. As a result, electric fields measured at the 
seafloor 6 by a receiving antenna 24 at a suitable horizontal range are dominated by 

1 0 the components of the source fields which have propagated downwards through the 
seafloor 6, along within the underlying strata 8, 12, and back up to the receiving 
antenna 24. Both the amplitude and the phase of the received signal depend on the 
resistivity structure of the underlying strata 8, 12 - and so, in principal, a survey 
consisting of many transmitter (source) and receiver (detector) locations can provide a 

15 multi-dimensional image, by geophysical inversion, of sub-seafloor resistivity. 

The technique described here exploits the large resistivity contrast that exists 
between a hydrocarbon reservoir (typically tens of Qm or higher) and the over- and 
under-lying sediments (typically ~ 2 Qm or less). Such a contrast has a detectable 
influence on controlled source electromagnetic data collected at the seafloor 6 above 

20 the reservoir 12. The effect of the reservoir is most detectable in controlled source 
electromagnetic data at an appropriate frequency, and if the horizontal range from 
source 22 to receiver 24 is of the order of 2 to 5 times the depth of burial of the 
reservoir 12 in typical situations. 

This following text describes specific geometric and data reduction approaches 

25 that allow the effect of a hydrocarbon reservoir on the outcome of a controlled source 
electromagnetic survey to be detected and analysed in practice. 
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Use of a mobile source and multiple fixed receivers 

In order to achieve a satisfactory survey outcome, it is essential to make 
controlled source electromagnetic measurements over a broad range of survey 
5 geometries - in other words, many transmitting locations and many receiving 
locations. The transmitter, i.e. the transmitting antenna, requires significant power to 
drive it, of the order tens of kilowatts or greater for signals detectable at ranges of 
several kilometres. It must therefore be connected by an umbilical cable 16 to the 
survey vessel 14 and this makes it relatively straightforward to make the transmitter 

10 mobile. It can then be towed in an appropriate survey pattern from the surface survey 
vessel 14. Since in many situations surveys of this kind are liable to take place over 
areas of the seafloor 6 where sensitive engineering installations exist or are planned, 
there are significant advantages to using a transmitter which does not physically come 
into contact with the seafloor 6. Provided that the separation between the transmitting 

15 antenna 22 and the seafloor 6 is small compared to a skin depth of the investigating 
field in seawater, the survey can still be completed satisfactorily. As a result, the most 
appropriate form of transmitter to use for this type of survey is one in which a 
neutrally buoyant horizontal electric dipole antenna 22 is towed behind the deep- 
towed vehicle 1 8 at a height of a few metres to a few tens of metres above the seafloor 

20 6. 

In the case of the receivers, there is an important advantage in using a static 
recording device. It is necessary to measure the alternating electric or magnetic field at 
the seafloor 6, resulting from the signal emitted by the transmitter. In typical 
applications, the signal-to-noise ratio of the received signal will be critical to the 
25 success and resolution of the survey, and so must be maximised. Moving the receiver 
inevitably generates noise, whether the signal measured is magnetic or electric field. 
In the case of electric fields, any motion of the receiver through the conducting 
seawater medium 4 in the presence of earth's geomagnetic field will generate an 
electromotive force across the receiving antenna 24. Receiver movements will 
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therefore map into spurious electric field signals in the recorded data. In the case of 
magnetic field recordings, there are also significant disadvantages to moving the 
receiver. Most importantly, if vector rather than scalar magnetometers are used (i.e. 
measuring individual directional components of the magnetic field), any variation in 

5 the orientation of the receiving antenna 24 will again lead to significant spurious 
signals, since the magnetic detecting element will detect changes in the component of 
the geomagnetic field aligned with it. As a consequence of these two effects, any 
translational movement of an electric field sensor or rotational movement of a 
magnetic field sensor will result in contamination of the received signal by motionally 

10 induced noise. 

For these reasons, it is desirable to carry out a controlled source 
electromagnetic survey to investigate or detect the presence of buried hydrocarbons 
using a combination of a mobile horizontal electric dipole source, equipped with a 
neutrally buoyant streamed antenna 22 and operated just above the seafloor 6; and an 

15 array of electric and/or magnetic field sensing receivers 24 placed statically on the 
seafloor 6 as shown in Figure 1. The receiving instruments 26 can be recovered using 
acoustically actuated release mechanisms to separate them from their ballast weights, 
allowing them to return to the sea surface for recovery, following standard 
oceanographic and marine geophysical practice. 

20 

The effects of propagation through the atmosphere 

There is an important additional factor that is crucial to the success or 
otherwise of the technique described here. In shallow water depths, it is possible for 
25 signals from the transmitter to follow a propagation path upwards through the water 
column to the surface; horizontally through the air; and back down through the water 
column to the seafloor receiver. This 'air wave' component contains no information 
about the sub-seafloor resistivity. It tends to dominate the received signal in shallow 
water and at long source-to-receiver offsets. The effect of the air wave can be 
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minimised by choosing appropriate transmission frequencies, and by targeting surveys 
on prospects in deep water and in which the target is at a relatively shallow depth 
below the seafloor. 

Figures 2 A, 2B and 2C show three cases of one-dimensional (ID) modelling 
5 showing signal amplitude as a function of source to receiver range and frequency. 
Each of the figures corresponds to a different depth of seawater overlying the seafloor. 

Figure 2A shows the modelled decay of electric field strength E as the distance 
R between the source (transmitter) and detector (receiver) increases. The detector is 
separated from the receiver along a line which runs parallel to, and passes through the 
10 transmitting dipole source. The model survey is undertaken above a semi-infinite 
subterranean strata of uniform resistivity 1 Qm, and in seawater of depth 750 m. The 
decay is shown for five different electromagnetic frequency components ranging from 
0.5 Hz to 8 Hz. 

Figure 2B represents the results of a similar modelling to that shown in Figure 

1 5 2 A, but for a seawater depth of 900 m. 

Figure 2C represents the corresponding results for a seawater depth of 1300 m. 
On all three sets of curves shown in Figures 2A, 2B and 2C, amplitude 
decreases rapidly with range, as expected. Additionally, amplitudes at a given range 
decrease with increasing frequency. This is because higher frequencies have shorter 

20 skin depths, and so experience increased attenuation. On each curve, the break in 
slope to a distinctly shallower gradient on the graph indicates the emergence of the air 
wave as the dominant signal at the receiver. For example, the curve corresponding to a 
signal component at 4 Hz shown in Figure 2B shows that the air wave becomes 
important in seawater of depth 900 m and using a 4 Hz source transmission frequency 

25 when the detector is more than about 3.5 km away from the source. It can be seen that 
this becomes an increasing problem in shallow water and at higher frequencies. For 
the technique described here to work most efficiently, the signal at the receiver must 
be primarily due to propagation through the seafloor, and not due to the air wave. For 
instance, for this condition to be met at 5 to 6 km source-detector offset for the model 
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subterranean strata described above, the frequency used for the survey should be no 
higher than 1 Hz in 1300 m water depth, or 0.5 Hz in 900 m water depth. 

Modelling can thus be used to plan a survey in order to determine the 
maximum distance allowable between source and detectors for various frequencies. 
5 Suitable detector deployment positions and source frequency ranges can thus be 
determined in advance of performing a survey. 

The use of survey result parameters derived from polarisation ellipse analysis 

10 Use of a horizontal electric dipole antenna produces electromagnetic fields at 

the seafloor that can be measured using electrometers or scalar or vector 
magnetometers. Scalar magnetometers are not used widely in practice. Current 
practice is to measure two or more orthogonal components of either electric or 
magnetic field. Up to six channels of data can usefully be recorded by the receiving 

15 instruments, corresponding to three orthogonal directional components each of 
electric and magnetic field. For operational reasons of instrument complexity and data 
storage, in most cases a sub-set of these is recorded. 

The commonest arrangement currently is to use an orthogonal pair of 
horizontal electric dipole receiving elements in each receiving instrument. In the 

20 horizontal plane, the transmitter generates both radial and azimuthal components of 
electric field at the receiver, and, since these have different amplitudes and phases in 
general, the electric field at the seafloor from a harmonic transmitter maps out a 
polarisation ellipse in the horizontal plane. The same applies in the case of horizontal 
magnetic field; and if vertical components are included, the planar ellipse becomes a 

25 polarisation ellipsoid. 

A straightforward approach to analysing controlled source electromagnetic 
data is to resolve the observed fields into radial and azimuthal components. However 
this suffers from contamination of field values by errors arising from small 
inaccuracies in source and receiver orientation and position. As an alternative, we 
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therefore propose the novel approach of analysing the data from the recorded 
components in terms of polarisation ellipse properties - e.g. the amplitude and/or 
phase of the component of the signal aligned along the major axis of the ellipse. Since 
this value is much less dependent on the accuracy of the source and receiver 

5 positioning - and in particular removes the significant errors that arise from even 
small uncertainties in receiver orientation - use of polarisation ellipse values can lead 
to smaller uncertainties in measured field properties at the receiver. We show below 
as an example that use of the amplitude along the major axis of the polarisation ellipse 
for horizontal electric field can be used as a robust measurement parameter for 

10 detecting the presence of buried hydrocarbon layers. In principle, other polarisation 
ellipse parameters can be used in a similar way. Parameters that could readily be 
measured include amplitude or phase along the major axis; the horizontal polarisation 
ellipse of either electric or magnetic field; or in either case a polarisation ellipsoid 
including a vertical component. 

15 

Presentation of survey results using normalised values 

As seen in Figures 2A, 2B and 2C above, a feature of a controlled source 
electromagnetic survey over a sedimented seafloor is that the received field amplitude 

20 drops very rapidly with increasing range. The received signal properties also depend 
on the orientations of the source and receiving dipoles. In a typical survey, the signal 
amplitude is likely to vary by several orders of magnitude over the useful set of 
source-detector (transmitter-receiver) offsets; while the phase lag of the received 
signal also increases steadily with increasing offset. The effects of buried structure 

25 beneath the survey show up as deviations in the signal from that which would be 
expected if the sub-surface consisted only of a homogeneous, isotropic half space. 

In order to visualise more clearly the effects of buried structure on controlled 
source electromagnetic signal properties, it is convenient to normalise the observed 
signals with respect to some reference model. A simple model to use would consist of 
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a water layer of finite thickness above the survey of the true depth; and a 
homogeneous isotropic half space representing the seafloor below the survey. If 
appropriate, and if the relevant a priori information is available, a more complex 
reference model may be used, although it is desirable to use always the simplest 
5 reference model that can broadly represent the large scale background properties of 
the seafloor. 

The normalisation is carried out for amplitude by dividing the observed field 
by that calculated for the reference model, using the same source-detector geometry as 
for the observed data. In the case of phase, the normalisation is carried out by 

10 subtracting the phase lag calculated for the reference model from that of the observed 
data. We show below that use of normalised values calculated in this way can 
dramatically enhance the visual presentation of controlled source electromagnetic data 
from a survey over a hydrocarbon reservoir. Our illustrations use normalised 
amplitudes, but normalised phases could equally be used. 

15 As an alternative, normalisation could be based on the survey data itself, for 

example using data collected adjacent to the target. 

The varying physical response of a survey depending on source orientation geometry 

20 Overall, the most crucial factor for the success or otherwise of the controlled 

source electromagnetic technique in practical applications related to hydrocarbon 
reservoirs is related to survey geometry. 

Figure 3 is a schematic plan view from above illustrating a suitable co-ordinate 
system for describing the relative placement of a controlled source electromagnetic 

25 source 30 and an electromagnetic detector 32. The position of the detector 32 with 
respect to the source 30 is most suitably described in polar co-ordinates, with the 
centre of the source 30 providing the origin of the polar co-ordinate system. The 
position of the detector 32 is defined by an azimuthal angle 9 and a distance R. The 
detector 32 and source 30 are considered to be co-planar. The angle 0 is measured 
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clockwise from a line 33 passing through, and running parallel to, the source axis, as 
indicated in Figure 3 by the line marked G=0°. A detector placed along the line 33, i.e. 
such that is has an azimuthal angle 9 of 0°, is referred to as being in an end-on 
position. A detector with an azimuthal angle 9 of 90°, such that it lies on the line 34 
5 marked 9=90° in Figure 3, is referred to as being in a broadside position. The electric 
field strength at the detector is considered in terms of a radial component E p and an 
orthogonal azimuthal component Ee as indicated in Figure 3. 

Depending on the relative placements and orientations of sources and 
receivers, the physics of propagation of the signal through the seafloor - and so the net 

10 effect on the properties of the received signal - can be more or less sensitive to 
different classes of sub seafloor structure. As a simple illustration of this in the 
context of hydrocarbon surveys, we consider a case in which the sub-seafloor structure 
can be represented by a simple stack of horizontal layers (a useful first approximation 
for many situations in sedimentary basins). 

15 Figure 4A shows in schematic vertical section an exemplary subterranean 

strata configuration. A section of seafloor 42 lies beneath 800 m of seawater 40. The 
strata beneath the seafloor 42 comprise a 1 km thick overburden layer 44, representing 
sediments above a hydrocarbon reservoir layer 46. This overburden layer 44 has low 
resistivity, primarily due to aqueous saturation of pore spaces. The middle layer 46 

20 corresponds to a 100 m thick hydrocarbon reservoir and has a resistivity perhaps 100 
times greater than the overburden 44. This is due to the presence of the non- 
conducting hydrocarbon occupying much of the pore spaces. Below the hydrocarbon 
reservoir layer 46 is a sub-reservoir region 48 which has low resistivity due to its 
similarity to the overburden layer 44 and extends downwards for an effectively 

25 infinite distance. 

In the case of an electromagnetic survey, we can consider the differing 
behaviours of electric currents generated by the source resolved along horizontal and 
vertical directions. The relationships between the electric currents flowing in two 
adjacent regions of space are determined by both galvanic (i.e. direct transfer of 
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charge) and inductive effects. Since charge is conserved, current leaving one volume 
of the sub-surface strata 44, 46, 48 and arriving in an adjacent volume along the 
direction of current flow are related to each other by a galvanic mechanism. On the 
other hand, if two volumes are close together but separated from each other along a 

5 direction orthogonal to current flow, then the coupling between the currents flowing in 
the two volumes will be primarily inductive. 

Applying this to our model of a sub-seafloor structure shown in Figure 4A, we 
can infer that the effect of the thin but resistive hydrocarbon reservoir layer 46 on the 
survey results will depend strongly on the direction of flow of the currents generated 

10 by the source. If the current at the base of the overburden layer 44 is dominantly 
horizontal, then coupling between this layer and the deeper layers 46, 48 will be 
dominated by inductive effects. Although little current will flow in the hydrocarbon 
reservoir layer 46, the inductive coupling between the overburden 44 and the 
underlying layer 48 will be only mildly affected. Currents in the two conducting layers 

1 5 44, 48 (overburden and underlying layers) will therefore be similar to the case without 
a hydrocarbon reservoir layer 46. The results of such a survey will therefore be only 
weakly affected by the presence of the hydrocarbon reservoir layer 46. 

In contrast, if our survey generates significant vertical components of electric 
current flow in the overburden layer 44, then galvanic effects along the current path 

20 will be strongly affected by the thin resistive hydrocarbon reservoir layer 46. The 
resistive hydrocarbon reservoir layer 46 will tend to block the current flow. As a 
result, currents flowing in the deeper layer 48 will be much reduced; and the overall 
pattern of current flow in the overburden layer 44 will therefore be significantly 
altered. In such a survey situation, the presence of the hydrocarbon will significantly 

25 influence the outcome. 

Hence it becomes of paramount importance in designing a practical survey 
approach for detecting buried hydrocarbon reservoirs to distinguish between source 
(transmitter) and detector (receiver) geometries in which coupling between layers is 
largely inductive between current sheets in a horizontal plane (in which case the 
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survey has little sensitivity to the presence of the reservoir); and those in which a 
significant component of vertical current flow occurs, in which case blocking of the 
passage of this current flow by the reservoir leads to a survey which is strongly 
sensitive to the presence of the reservoir. 

5 

A specific approach to combining two contrasting source geometries when 
prospecting for hydrocarbon filled reservoirs 

By selecting a suitable survey geometry, it is possible to exploit the diverse 

10 properties of electromagnetic induction outlined above, by designing a survey in such 
a way that it provides data that are sensitive to the presence of thin resistive layers 
(exploiting largely galvanic effects and vertical components of induced current flow); 
while simultaneously obtaining data that are more sensitive to the much larger scale 
structure, but relatively insensitive to a thin hydrocarbon layer (exploiting the 

15 dominantly inductive effects). The reason why the latter component is essential is that 
many features other than hydrocarbon reservoirs can affect the resistivity beneath the 
seafloor and the results of a survey. For example, resistivity often increases steadily 
with depth in submarine sedimentary basins, due to the progressive expulsion of pore 
fluids by rising overburden pressure. 

20 Figure 4B shows in vertical section a highly schematic exemplary subterranean 

strata configuration which exhibits increasing resistivity with depth. A section of 
seafloor 52 lies beneath 800 m of seawater 50. The strata beneath the seafloor 52 
comprise a series of layers of sediment of increasing resistivity. A first layer 54 has a 
uniform resistivity of 1 Qm and a thickness of 1 km. A second layer 56 has a uniform 

25 resistivity of 5 Qm and a thickness of 1 km. A third layer 58 has a uniform resistivity 
of 10 Clm and a thickness of 1 km. A fourth layer 60 has a uniform resistivity of 50 
Qm and a thickness of 1 km. Beneath the fourth layer 60 is a fifth layer 62 which has 
a resistivity of 100 Qm and extends downwards for an infinite extent. 
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The increasing resistivity indicated in Figure 4B leads to longer skin depths, 
and hence to larger observed amplitudes and smaller observed phase lags for all 
detector placements. Although this does not exactly reproduce the same galvanic 
effects of a thin resistive layer, the overall effect on observed fields in the end-on 

5 geometry is likely to be very similar to (and difficult or impossible to distinguish 
from) the effect of a hydrocarbon reservoir. Designing a survey in such a way that 
different parts of the resulting data are more or less sensitive to the presence of certain 
key features is essential for removing ambiguity in the interpretation. 

Recent studies of volcanic systems at mid-ocean ridges and of sediments 

10 beneath resistive basalt layers (MacGregor et aL, 1998, 2001; MacGregor & Sinha, 
2000) have demonstrated that improved model resolution and reduced interpretational 
ambiguity can be achieved by using received electric field components from two 
distinct source-detector geometries. This approach is now applied to the case of 
surveying for thin hydrocarbon reservoirs. 

15 Figures 4C, 4D, 4E and 4F show survey simulation results for the two 

subterranean strata configurations shown in Figures 4A and 4B. The method relies on 
collecting survey data using two distinct geometric modes. The first mode 
corresponds to an end-on geometry, in which, as described above, the receiver lies 
along the axis of the transmitting dipole. The source-detector azimuth 0 as defined in 

20 Figure 3 is 0°, and the field observed at the receiver is dominated by the radial 
component E p . The second mode corresponds to a broadside geometry, in which, as 
described above, the source-detector azimuth 0 is 90°, and the field observed at the 
receiver is dominated by the azimuthal component Eq. 

For the case of hydrocarbon exploration, where the target is a thin resistive 

25 layer embedded between a more conductive underburden and overburden, such as that 
shown schematically in Figure 4A, the received signal depends on the azimuthal angle 
8 of the detector. Because of the transmitter field pattern, the end-on, or E p , geometry 
results in a significant contribution to the observed field at the seafloor by the vertical 
component of current flow. The broadside, or E e , geometry results in fields at the 
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seafloor that are more dependent on the contribution of inductively coupled currents 
flowing in horizontal planes. As a result, the presence of a resistive hydrocarbon 
reservoir layer produces a significant increase at certain ranges in the end-on (E p ) 
fields, while having virtually no effect on the broadside (E e ) fields. The result is 
5 'splitting' of amplitudes between the two modes, and this splitting is characteristic of 
a structure in which resistivity first increases, in this case due to a hydrocarbon layer, 
and then decreases with depth. 

Figure 4C is a graph showing the modelled amplitudes of electric field E along 
the major axes of the horizontal polarisation ellipses for two source-detector 

10 geometries and as a function of distance R between source and detector. These are 
determined for the model of sub-seafloor structure indicated in Figure 4A. Four curves 
are shown. The curve labelled E 0 shows the field strength as a function of distance for 
end-on geometries, the curve labelled E 9 o shows the same for broadside geometries. 
The curves labelled No and N 90 are the corresponding normalisation curves. These are 

15 determined for the same geometries as for Eo and E 90 , but represent the modelled 
amplitudes of electric field along the major axes of the horizontal polarisation ellipses 
where the subterranean strata comprise a homogeneous isotropic half space of 
resistivity 1 Qm, such as previously described. 

Figure 4D is a graph showing the set of curves which correspond to those 

20 curves shown in Figure 4C, but for the subterranean strata configuration shown in 
Figure 4B. 

Figure 4E is a graph showing the normalised values of the fields Eo and E90 
shown in Figure 4C. These fields have been normalised by the fields N 0 and N 90> also 
shown in Figure 4C, as discussed above. 
25 Figure 4F is a graph showing the normalised values of the fields Eo and E 9 o 

shown in Figure 4D. These fields have been normalised by the fields N 0 and N 90 , also 
shown in Figure 4D, as discussed above. 

It can be seen in Figure 4E that the presence of the reservoir has relatively little 
effect on the broadside amplitudes (i.e. the curve labelled E 90 /N 90 is close to unity), 
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but that between 4 and 6 km offset, it produces a substantial increase in the end-on 
amplitudes (i.e. the curve labelled Eo/N 0 is strongly enhanced over this range). The 
hydrocarbon layer in this particular model leads to amplitude splitting between the 
two geometrical modes by about a factor of 10. 

5 It can be seen in Figure 4F that the steadily increasing resistivity with depth 

model shown in Figure 4B strongly affects both geometries (i.e. both curves are 
strongly enhanced between 4 and 6 km offset) and no significant splitting between the 
two geometric modes occurs. 

The end-on geometry data alone cannot distinguish between the two 

10 subterranean strata configurations shown in Figures 4A and 4B. These comparative 
examples conclusively demonstrate the fact that collection of survey data from both 
end-on and broadside geometries is needed, in order to reliably detect the hydrocarbon 
layer and eliminate "false-positive" detection of a rock formation containing no 
hydrocarbon reservoir, but merely having an increasing resistivity with depth, which is 

15 not unusual. 

Lastly, the normalised curves in Figures 4E and 4F start off close to unity, at 
short ranges where the data are insensitive to the buried structure. The effects of the 
buried structure become greatest at between 5 and 6 km range. Beyond this, the 
normalised curves rapidly return towards a value of 1, because at long ranges, air 
20 wave propagation (insensitive to sub-seafloor structure) begins to dominate, thus 
masking the effects of the sub surface. As a result, the presence of the target structure 
can best be observed in controlled source electromagnetic data of this type over a 
relatively narrow window of source-detector separation ranges. 

An efficient electromagnetic survey may use a large number of detectors 
25 distributed over a target area and a signal transmitter which operates continuously as it 
is towed along an extended tow path. This follows from the relatively low cost of 
deploying detectors and the relatively high cost of deploying a transmitter. 

Figure 5 is a schematic plan view showing an example layout of sixteen 
detectors 64 distributed across a section of seafloor 65 above subterranean reservoir 
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66. The reservoir 66 has a linear extent of several km and its boundary is indicated by 
a heavy line 67. In this example, the detectors 64 are uniformly distributed in a 
square-grid pattern so as to approximately cover the substantially square reservoir 66. 
In performing a survey, a source (not shown) starts from the location marked 'A' in 
5 Figure 5 and is towed, whilst broadcasting continuously, along a path indicated by the 
broken line 68, the survey is completed when the source reaches the location marked 
*B\ Data are continuously recorded by the detectors 64 throughout the towing process 
and the position of the source transmitter relative to the detector network is also 
logged. 

10 During the towing process, each of the detectors 64 presents several different 

orientations to the source. For example, when the source is directly above the detector 
position Dl and on the vertical section of the tow path, the detectors at positions D2 
and D3 are at different ranges in an end-on position, the detectors at positions D4 and 
D5 are at different ranges in a broadside position and the detector at position D6 is 

1 5 midway between. However, when the source later passes over the detector position Dl 
when on the horizontal section of the tow path, the detectors at positions D2 and D3 
are now in a broadside position, and the detectors at position D4 and D5 are in an end- 
on position. Thus, in the course of a survey, and in conjunction with the positional 
information of the source, data from the detectors 64 can be used to provide details of 

20 the signal transmission through the subterranean strata for a comprehensive range of 
distances and orientations between source and detector, each with varying galvanic 
and inductive contributions to the signal propagation. In this way a simple continuous 
towing of the source transmitter can provide a detailed survey which covers the extent 
of the subterranean reservoir 66. 

25 Although the above example is based on a square detector grid, it will be 

understood that a wide variety of detector placements may be used, for example other 
high symmetry regular grids, such as triangular or rectangular, may be used. In 
addition irregular grids may be used that have no high level of symmetry. 
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On the detection of the edges of a reservoir b y this technique 

We showed above that our proposed arrangement of survey geometries can 
both detect a hydrocarbon reservoir, and distinguish it from a general increase in 
5 resistivity with depth. The modelling though was based only on 1-D assumptions, i.e. 
horizontal layers of infinite horizontal extent. By using higher dimensional modelling, 
we show below that the survey technique also works well in the case of reservoirs of 
finite extent; and indeed that it can be used as a reliable means of locating the edges of 
a reservoir structure. 

10 Figure 6 A shows a schematic section view of a model subterranean strata 

configuration containing a hydrocarbon reservoir layer of finite horizontal extent. A 
section of seafloor 70 lies beneath a body of sea 72 which is 1.15 km deep and has a 
resistivity of 0.31 Qm. A planar hydrocarbon reservoir 74 is 0.15 km thick and its 
base is 1 km below the seafloor 70. The reservoir 74 has a resistivity of 100 Qm and 

15 extends infinitely out of the plane of Figure 6A, and semi-infinitely within the plane 
of Figure 6A. A uniform sedimentary structure 76 of resistivity 1 Clm surrounds the 
reservoir 74 and fills the remaining space beneath the seafloor 70. For the purposes of 
the model, two transmitter positions are considered. A first dipole transmitter 75 is 
aligned parallel to and directly above the edge 77 of the reservoir 74 on the seafloor 

20 70. The dipole axis of the transmitter 75 is perpendicular to the plane of Figure 6A 
and a cross is used to simply indicate its position rather than represent its structure. A 
second transmitter 78 is also on the seafloor 70 and aligned parallel to the first 
transmitter 75. The second transmitter 78 is horizontally displaced from the first 
transmitter 75 by 1 km in a direction which positions it directly above the reservoir 74 

25 as indicated in Figure 6A. 

Figures 6B and 6C show the results of 2.5 D modelling (3-D source, 2-D 
resistivity structure) of the exemplary subterranean strata configuration shown in 
Figure 6A for the first and second transmitter positions 75 and 78 respectively in order 
to investigate the effects of the finite areal extent of the resistive hydrocarbon 
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reservoir 74. Figures 6B and 6C are plan views representing the normalised fields at 
the seafloor 70 as a function of position. The normalised field component shown in 
each case is the semi-major axis of the horizontal polarisation ellipse of electric field. 
In Figure 6B the first transmitter 75 is at the origin of the local co-ordinate system and 
5 its dipole axis is parallel to the vertical axis in this system. In Figure 6C the second 
transmitter 78 is at the origin of the local co-ordinate system and its dipole axis is 
parallel to the vertical axis in this system. 

The normalised field amplitude is shown in a grey scale representation, while 
contour lines indicate the absolute value of amplitude. Figures 6B and 6C show that 
10 the reservoir structure causes substantial increases in the fields at positions close to 
end-on geometries, this is especially so when the source and detector are both over the 
reservoir, as in Figure 6C. If the source or detector are at or outside the edge of the 
reservoir, the amplitude increase is not seen. At detector locations close to broadside 
geometry, the increase in amplitude is again not seen, even if both source and detector 
15 are over the reservoir. Figures 6B and 6C illustrate that the amplitude splitting effect 
is observable over reservoirs of finite size, as well as over the 1-D (layered) structures 
discussed previously. 

The marked difference between Figures 6B and 6C shows that the observed 
fields are extremely sensitive to the location of the source (transmitter) with respect to 
20 the edge of the reservoir. In Figure 6B, the marked asymmetry of the fields about the y 
axis shows that - provided the source is over the reservoir - the observed fields are 
also highly sensitive to the location of the detector (receiver) with respect to the edge 
of the reservoir. Thus by careful siting of both source tow tracks and detectors relative 
to a suspected hydrocarbon bearing structure, the survey method described here can 
25 provide detailed information on the areal extent of the hydrocarbon, as well as on its 
existence or otherwise. 

In addition to initially surveying for subterranean hydrocarbon reservoirs, the 
utility of the invention in detecting the edges of reservoirs makes it applicable to 
assessing changes in hydrocarbon reservoir content over a period of time. This is of 
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particular value, for example, where a hydrocarbon reservoir is being exploited. As 
water replaces the hydrocarbon which is drawn from the reservoir, regular follow up 
surveys can be used to determine the change in the distribution of the remaining 
hydrocarbon. Such surveys could therefore be used to characterise the evolution of 
5 the reservoir through time, and hence enable more efficient extraction and reservoir 
management. 

Source and receiver locations - design considerations for real surveys 

10 The methodology described in this document can be readily applied to real 

surveys for hydrocarbon resources on continental margins. The required transmitter 
and receiver characteristics in terms of power, signal to noise ratios and operating 
parameters can be met by existing instrumental technology. 

In order to plan a successful survey, the following factors are important: 

15 An understanding of the need to collect observations over an appropriate set of 

source-detector separation ranges, and using both broadside and end-on surveying 
geometries; 

A designed layout of the detector array and of source tow line locations and 
orientations optimised for the above; 
20 A designed layout for the survey that also takes account of the need to locate 

the limits of the areal extent of any reservoir, based on the approach presented in 
Figures 6A, 6B and 6C; 

Prior modelling to establish the range of transmission frequencies, and the set 
of source-detector separation ranges, that should be used, taking into consideration the 
25 expected resistivity structure of the seafloor, the depth of the target, the depth of the 
seafloor, and the influence of the air wave on the data collected. 
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Conclusions 

We have shown that the application of a specific set of methods to a controlled 
source electromagnetic survey of the seafloor can allow the technique to be 
5 successfully applied to the problem of detecting hydrocarbons beneath the seafloor. 
Specifically, we have shown that: 

The survey should be carried out using a combination of a mobile horizontal 
electric dipole source transmitter, equipped with a neutrally buoyant streamed antenna 
and operated just above the seafloor; and an array of electric and/or magnetic field 
10 sensing receivers placed statically on the seafloor. 

The effects of propagation of the signal through the atmosphere (the e air 
wave') are significant at high frequencies, at long source-detector separation 
distances, and in shallow water. This effect limits the applicability of the method by 
favouring detection of structures relatively shallowly situated beneath the seafloor, but 
15 in deep water; and limits the choice of frequencies and source-detector offsets for the 
survey. 

Improvements in data analysis can be achieved by making use of field 
measurement parameters aligned along the major axis of the polarisation ellipse at the 
seafloor. 

20 For purposes of data presentation and interpretation it is desirable to make use 

of field parameters that have been normalised with reference to an appropriate 

simplified model of the sub-seafloor structure. 

In order to resolve the presence of any hydrocarbon, and to distinguish its 

effect on the data from other likely structures and reduce the ambiguity of 
25 interpretation, it is essential to collect survey data from both the end-on and broadside 

geometric configurations; and to analyse the data in terms of splitting between the two 

geometric modes. 

2.5 D modelling shows that provided an appropriately designed array of 
receivers, and an appropriate set of transmitter tow tracks, are used, then the method 
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can yield valuable information about the limits of the areal extent of any sub-seafloor 
hydrocarbon reservoir, as well as detecting its presence. 

Finally it will be understood that the invention is equally applicable to 
surveying of freshwater, for example large lakes, so that references to seafloor, 
5 seawater etc. should not be regarded as limiting. 
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CLAIMS 

1. An electromagnetic survey method for surveying an area previously identified 
as potentially containing a subsea hydrocarbon reservoir, comprising: 

5 obtaining first and second survey data sets with an electromagnetic source 

aligned end-on relative to a first detector and aligned broadside relative to a second 
detector. 

2. An electromagnetic survey method according to claim 1, wherein the first and 
10 second survey data sets are obtained concurrently. 

3. An electromagnetic survey method according to claim 1 or claim 2, wherein a 
right angle is formed between a first line leading from the first detector to the 
electromagnetic source and a second line leading from the second detector to the 

15 electromagnetic source. 

4. An electromagnetic survey method for surveying an area previously identified 
as potentially containing a subsea hydrocarbon reservoir, comprising: 

obtaining first and second survey data sets with an electromagnetic source 
20 aligned end-on and broadside respectively relative to a first detector. 

5. An electromagnetic survey method according to claim 4, wherein the first and 
second survey data sets are obtained consecutively. 

25 6. An electromagnetic survey method according to claim 4 or claim 5, wherein 
the electromagnetic source is moved in a direction along a line connecting the 
electromagnetic source to the detector to obtain one of the first and second survey data 
sets and again transverse thereto to obtain the other of the first and second survey data 
sets. 
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7. An electromagnetic survey method according to any one of the preceding 
claims, wherein the first and second survey data sets are obtained by moving the 
electromagnetic source relative to each detector. 

5 

8. An electromagnetic survey method according to claim 7, wherein the 
electromagnetic source is a mobile horizontal electric dipole source equipped with a 
streamed antenna which is towed above the seafloor to obtain the first and second 
survey data sets. 

10 

9. An electromagnetic survey method according to any one of the preceding 
claims, wherein each detector is static while the first and second survey data sets are 
obtained. 

15 10. An electromagnetic survey method according to claim 9, wherein each 
detector is placed on the seafloor while the first and second survey data sets are 
obtained. 

11. A method of analysing results from an electromagnetic survey of an area 
20 potentially containing a subsea hydrocarbon reservoir, comprising: 

providing first and second survey data sets obtained from an electromagnetic 
source aligned respectively end-on and broadside relative to a detector; and 

combining the first and second survey data sets to obtain a results data set that 
represents a difference between the end-on and broadside alignments. 

12. A method of analysing results from an electromagnetic survey according to 
claim 11, further comprising: 
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normalising each of the first and second survey data sets relative to respective 
first and second normalisation data sets or functions specific to the end-on and 
broadside alignments respectively, prior to the combining. 

5 13. A method of analysing results from an electromagnetic survey according to 
claim 12, wherein the first and second normalisation data sets or functions are 
calculated from a rock formation model. 

14. A method of analysing results from an electromagnetic survey according to 
10 claim 12, wherein the first and second normalisation data sets or functions are 

calculated from the first and second survey data sets. 

15. A method of analysing results from an electromagnetic survey according to 
any one of claims 1 1 to 14, wherein the first and second data sets each comprise radial 

15 and azimuthal components of electric field or magnetic field measured at the detector, 
the method further comprising: 

transforming the radial and azimuthal components into at least one 
polarisation ellipse parameter, prior to the combining. 

20 16. A method of analysing results from an electromagnetic survey according to 
claim 15, wherein the at least one polarisation ellipse parameter is the amplitude 
and/or phase of the component of the electric field or magnetic field aligned along a 
major axis of the ellipse. 

25 17. A method of analysing results from an electromagnetic survey according to 
any one of claims 1 1 to 1 6, further comprising: 

visually representing the results data set in a plot of at least two dimensions 
corresponding to the survey area. 
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18. A method of analysing results from an electromagnetic survey according to 
claim 17, wherein the plot includes markings of areas of equal or similar 
electromagnetic field strength. 

5 19. A method of analysing results from an electromagnetic survey according to 
claim 18, when appended to any one of claims 12 to 14, wherein the plot includes 
lines of equal absolute electromagnetic field strength. 

20. A computer program product bearing machine readable instructions for 
10 implementing a method of analysing results from an electromagnetic survey according 

to any one of claims 1 1 to 19. 

21. A computer apparatus loaded with machine readable instructions for 
implementing the method of analysing results from an electromagnetic survey 

15 according to any one of claims 11 to 19. 

22. A method of planning an electromagnetic survey of an area identified as 
potentially containing a subsea hydrocarbon reservoir, comprising: 

creating a model of the area to be surveyed, including a rock formation 
20 containing a hydrocarbon reservoir and a body of water above the rock formation; 

setting values for water depth, depth below the seafloor of the hydrocarbon 
reservoir, and resistivity structure of the rock formation; 

performing a simulation of an electromagnetic survey in the model of the 
survey area by calculating first and second survey data sets for an electromagnetic 
25 source aligned end-on and broadside relative to a detector; and 

combining the first and second survey data sets to obtain a results data set that 
represents a difference between the end-on and broadside alignments. 
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23. A method of planning an electromagnetic survey according to claim 22, further 
comprising: 

repeating the simulation for a number of source frequencies and/or source-to- 
detector distances in order to select optimum surveying conditions in terms of source 
5 frequency and source-to-detector distance for probing the hydrocarbon reservoir. 

24. A method of planning an electromagnetic survey according to claim 22 or 23, 
wherein the model includes a body of air above the body of water, and wherein the 
simulation takes account of signal propagation paths including the body of air when 

1 0 calculating the first and second survey data sets. 

25. A method of planning an electromagnetic survey according to any one of 
claims 22 to 24, further comprising: 

normalising each of the first and second survey data sets relative to respective 
15 first and second normalisation data sets or functions specific to the end-on and 
broadside alignments respectively, prior to the combining. 

26. A method of planning an electromagnetic survey according to any one of 
claims 22 to 25, wherein the first and second data sets each comprise radial and 

20 azimuthal components of electric field or magnetic field measured at the detector, the 
method further comprising: 

transforming the radial and azimuthal components into at least one 
polarisation ellipse parameter, prior to the combining. 

25 27. A method of planning an electromagnetic survey according to any one of 
claims 22 to 26, further comprising visually representing the results data set in a plot 
of at least two dimensions corresponding to the survey area. 
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10 



28. A method of planning an electromagnetic survey according to claim 27, 
wherein the plot includes markings of areas of equal or similar electromagnetic field. 

29. A method of planning an electromagnetic survey according to claim 28, when 
appended to claim 25, wherein the plot includes lines of equal absolute 
electromagnetic field strength. 

30. A computer program product bearing machine readable instructions for 
implementing the method of planning an electromagnetic survey according to any one 
of claims 22 to 29. 

31. A computer apparatus loaded with machine readable instructions for 
implementing the method of planning an electromagnetic survey according to any one 
of claims 22 to 29. 
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